AD-A212  202 


nri"  r  -  ,  r;y 

U  I  l  ;.}  .  .1-1 

I  David  Taylor  Research  Center 

Bethesda,  Maryland  20084-5000 


DTRC/CMLD-S9/07  July  1989 

Computation,  Mathematics  and  Logistics  Department 
Research  and  Development  Report 


A  Mechanical  Power  Flow  Capability 
for  the  Finite  Element  Code  NASTRAN 

by 

Stephen  A.  Hambric 


Approved  for  public  release;  distribution  unlimited. 


89  9  13  103 


CODE  Oil  DIRECTOR  OF  TECHNOLOGY,  PUNS  AND  ASSESSMENT 
12  SHIP  SYSTEMS  INTEGRATION  DEPARTMENT 

14  SHIP  ELECTROMAGNETIC  SIGNATURES  DEPARTMENT 

15  SHIP  HYDROMECHANICS  DEPARTMENT 

16  AVIATION  DEPARTMENT 

17  SHIP  STRUCTURES  AND  PROTECTION  DEPARTMENT 

18  COMPUTATION,  MATHEMATICS  &  LOGISTICS  DEPARTMENT 

19  SHIP  ACOUSTICS  DEPARTMENT 

27  PROPULSION  AND  AUXILIARY  SYSTEMS  DEPARTMENT 

28  SHIP  MATERIALS  ENGINEERING  DEPARTMENT 


DTRC  ISSUES  THREE  TYPES  OF  REPORTS: 

_  * 

1.  DTRC  reports,  a  formal  series,  contain  information  of  permanent  technical  value. 
They  carry  a  consecutive  numerical  identification  regardless  of  their  classification  or  the 
originating  department. 

2.  Departmental  reports,  a  semiformal  series,  contain  information  of  a  preliminary, 
temporary,  or  proprietary  nature  or  of  limited  interest  or  significance.  They  carry  a 
departmental  alphanumerical  identification. 

3.  Technical  memoranda,  an  informal  series,  contain  technical  documentation  of 
limited  use  and  interest.  They  are  primarily  working  papers  intended  for  internal  use.  They 
carry  an  identifying  number  which  indicates  their  type  and  the  numerical  code  of  the 
originating  department.  Any  distribution  outside  DTRC  must  be  approved  by  the  head  of 
the  originating  department  on  a  case-by-case  basis. 


UNCLASSIFIED 


U  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2*  SECURITY  CLASSIFICATION  AUTHORITY 


REPORT  DOCUMENTATION  PAGE 


lb  RESTRICTIVE  MARKINGS 


2b  DECLASSIFICATION  / DOWNGRADING  SCHEDULE 


3  DISTRIBUTION /AVAILABILITY  Of  REPORT 

Approved  for  public  release;  distr ibut ion 

unlimited 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

CMLD-89/07 


S  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6«  NAME  OF  PERFORMING  ORGANIZATION  6b  OFFICE  SYMBOL 

David  Taylor  Research  Center  W eppbcsble) 
_  Code  1844 


6c  ADDRESS  (City.  Stete.  end  HP  Code) 

Bethesda,  MD  20084-5000 


7«.  NAME  OF  MONITORING  ORGANIZATION 


7b.  AOORESS  (City,  S»4»*.  end  ZIP  Code) 


84  NAME  OF  FUNOING/ SPONSORING 
ORGANIZATION 

Office  of  Naval  Technology 


8c  AOORESS  (City.  St4t*.  end  ZIP  Cod*) 

800  N.  Quincy  Street 
Arlington,  VA  22217 


ll  fir t£  (Include  Security  Clesuhcetion) 


8b.  OFFICE  SYMBOL 
Of  epplkeble)  ■ 

0NT 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


10  SOURCE  OF  FUNOING  NUMBERS 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO 

NO 

NO 

ACCESSION  NO 

62323N 

R2320 

l- 1903-00  l 

12  PERSONAL  AUTHOR(S) 

_ Stonhen  A.  Haabc 

134  Type  of  REPORT 


|Ub  time  COVEREO 
I  FROM _  TO 


|14  DATE  OF  REPORT  (Y»4C.  Month.  D*y)  |1 5  PAGE  COUNT 

I  1989,  Jul v  12 


COSATl  COOES 


GROUP  I  SUB-GROUP 


18.  SUBJECT  TERMS  (Confmu*  on  reverse  if  iwniuiy  4 nd  identify  by  block  number) 

>  Power  flow  Mechanical  Intensity 
Finite  Element  Method 
Structural  Dynamics  . 


19  ABSTRACT  (Continue  on  reverse  if  necessity  end  identity  by  block  number) 

The  identification  of  power  flow  paths  in  dynamically  loaded  structures  is  an  important,  but 
currently  unavailable,  capability  for  the  finite  element  analyst.  For  this  reason,  methods  for 
calculating  power  flows  and  mechanical  intensities  in  finite  element  models  are  developed  here. 
Formulations  for  calculating  input  and  output  powers,  power  flows,  mechanical  intensities,  and 
power  dissipations  for  beam,  plate,  and  solid  element  types  are  -der-ived .  NASTRAN  is  used  to 
calculate  the  required  velocity,  force,  and  stress  results  of  an  analysis,  which  a  post-proces¬ 
sor  then  uses  to  calculate  power  flow  quantities.  The  SDRC  1-deas  Super tatr  module  is  used  to 
view  the  final  results.  Test  models  include  a  simple  truss  and  a  beam-stiffened  cantilever 
plate.  Both  test  cases  showed  reasonable  power  flow  fields  over  low  to  medium  frequencies. 
Future  work  will  include  testing  with  more  complex  models,  developing  an  interactive  graphics 
program  to  view  easily  and  efficiently  the  analysis  results,  applying  shape  optimization 
methods  to  the  problem  with  power  flow  variables  as  design  constraints,  and  adding  tiie  power 
flow  capability  to  NASTRAN. 


21  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


zo  distribution /availability  of  abstract 
Qunclassifieo/unlimited  □  same  as  rpt  □  otic  users 


22b  telephone  (include  Are*  Code)  22 1  OFFICE  SYMBOL 

(301)227-1645 


00  fORM  1473,84 MAR  B 3  APR tbition  W4y  b* um4 until •■RRuttad  SECURITY  CLASSIFICATION  Of  This 

All  otbei  editions  ere  obsolete  '  UNCLASSIFIED 


SECURITY  CLASSIFICATION  Of  This  page 


CONTENTS 


Page 


ABSTRACT .  1 

ADMINISTRATIVE  INFORMATION .  1 

INTRODUCTION .  1 

THE  FINITE  ELEMENT  SOLUTION .  3 

General  Methods .  3 

Power  Flow  and  Mechanical  Intensity .  4 

Damping  and  Power  Dissipation .  5 

Power  Input .  7 

Power  Output .  7 

Power  Balance  . . 7 


Element  Formulations  . 

Beam  Elements . 

Power  Flow  and  Mechanical  Intensity 

Power  Dissipation . 

Plate  Elements . 

Power  Flow  and  Mechanical  Intensity 

Power  Dissipation . 

Solid  Elements . 

Power  Flow  and  Mechanical  Intensity 

Power  Dissipation . 

Scalar  Elements . 


8 

8 

8 

10 

10 

10 

12 

12 

12 

13 

13 


Computer  Methods  . .  13 

NASTRAN  Solution .  14 

Power  Flow  Algorithm .  14 

Post-Processing .  15 

TEST  CASES .  15 


Simple  Truss .  15 

Problem  Statement .  15 

Results  .  16 

Beam-Stiffened  Cantilever  Plate .  21 

Problem  Statement .  21 

Results  .  22 


SUMMARY  AND  FUTURE  WORK .  25 

APPENDIX:  NASTRAN  INPUT  DATA .  28 

REFERENCES .  29 


iii 


FIGURES 


Fig.  1.  Sample  Power  Flow  Diagram . 

Fig.  2.  The  BAR  Element . 

Fig.  3.  The  QUAD  Element . 

Fig.  4.  Power  Flow  Solution  Process . 

Fig.  5.  Simple  Truss  Problem . 

Fig.  6.  Power  Flows  and  Dissipations  for  a  Single  Frequency  .  . 

Fig.  7.  Spectrum  Plots  of  Power  Flows  for  Truss  Section  1 .  .  .  . 

Fig.  8.  Spectrum  Plots  of  Power  Flows  for  Truss  Sections  2  and 
3  . 

Fig.  9.  Power  Flow  Diagrams  for  Truss  Problem . 

Fig.  10.  Beam-Stiffened  Cantilever  Plate  Problem . 

Fig.  11.  Input  and  Output  Powers  for  Cantilever  Plate . 

Fig.  12.  Power  Flows  for  Three  Locations  Along  Beam 

Stiffener . 

Fig.  13.  Mechanical  Intensity  Magnitudes,  f=7.55  Hz . 

Fig.  14.  Mechanical  Intensity  Directions  £=7.55  Hz . 


4 

9 
11 
14 
16 
17 
]  s 

19 

20 


24 

25 

26 


j  *>•.  ■  -;.'i  ;v  '('os 


i  M 

"  .1 

i 


iv 


ABSTRACT 


The  identification  of  power  flow  paths  in  dynamically 
loaded  structures  is  an  important ,  but  currently  unavailable, 
capability  for  the  finite  element  analyst.  For  this  reason, 
methods  for  calculating  power  flows  and  mechanical 
intensities  in  finite  element  models  are  developed  here. 
Formulations  for  calculating  input  and  output  powers,  power 
flows,  mechanical  intensities,  and  power  dissipations  for 
beam,  plate,  and  solid  element  types  are  derived.  NASTRAN 
is  used  to  calculate  the  required  velocity,  force,  and  stress 
results  of  an  analysis,  which  a  post-processor  then  uses  to 
calculate  power  flow  quantities.  The  SDRC  I-deas  Supertab 
module  is  used  to  view  the  final  results.  Test  models  include 
a  simple  truss  and  a  beam-stiffened  cantilever  plate.  Both  test 
cases  showed  reasonable  power  flow  fields  over  low  to 
medium  frequencies.  Future  work  will  include  testing  with 
more  complex  models,  developing  an  interactive  graphics 
program  to  view  easily  and  efficiently  the  analysis  results, 
applying  shape  optimization  methods  to  the  problem  with 
power  flow  variables  as  design  constraints,  and  adding  the 
power  flow  capability  to  NASTRAN. 

ADMINISTRATIVE  INFORMATION 

This  work  was  performed  under  Program  Element  62323N,  Task  Area 
R2320,  Acoustic  Silencing  Task  (Auxiliary  Systems),  and  Work  Unit  1-1903- 
001-10.  The  DTRC  project  managers  were  D.E.  Goldsmith  and  L.T.  Mo. 
DTRC  Code  274. 


INTRODUCTION 

Structure-borne  sound  is  the  vibrational  energy  which  travels  through 
dynamically  loaded  mechanical  systems.  This  vibrational  energy  is  radiated 
eventually  into  an  acoustic  medium  as  noise.  An  example  cited  by  Wohlever 
and  Bernhard1  is  an  airplane  wing  loaded  by  engine  vibrations.  The  vibrational 
energy  travels  along  the  wing  to  the  fuselage  and  is  radiated  as  sound  into  the 
cabin.  Architects  face  the  problem  of  structure-borne  sound  in  hotels  and 
apartment  buildings,  where  vibrational  energy  flows  through  walls  and  floors 
and  is  radiated  as  sound  into  other  rooms.  This  problem  is  addressed  by 
Luzzato  and  Or  to  la. 2 

Dynamic  analyses  performed  to  solve  these  problems  produce  large 
amounts  of  data.  The  analyst  is  then  faced  with  the  problem  of  interpreting 
the  output.  Tabular  printouts,  spectrum  plots,  and  deformed  shape  plots  are 
ail  useful  ways  of  defining  the  state  of  a  structure.  Another  way  of  quantifying 
the  propagation  of  structure-borne  sound  is  the  calculation  of  power  Hows. 
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This  method  identifies  the  magnitude  and  direction  ot'  the  power  at  any 
location  in  a  structure  and  helps  an  analyst  find  the  dominant  paths  of  energy 
flow  and  the  energy  sinks  in  a  given  problem.  Understanding  the  paths  of 
energy  which  flow  from  a  vibration  source  to  other  parts  of  a  structure  helps 
an  engineer  to  pinpoint  and  correct  vibration  problems. 

For  naval  applications,  such  as  machinery  in  submarines  and  surface 
ships,  quantifying  power  flow  paths  could  lead  to  reduction  of  noise.  Applying 
isolators  and/or  active  control  devices  to  a  machinery  system  could  reduce 
power  flows  along  particular  paths  (to  the  hull  for  example)  and  therefore  also 
reduce  sound  radiation. 

The  following  important  terms  are  used  in  this  study:  power  flow, 
which  is  actually  power  or  energy  flow  but  is  termed  power  flow  by  the 
scientific  community  of  this  field;  mechanical  intensity,  which  is  power  flow 
per  unit  area;  and  power  dissipation,  which  is  the  time  rate  of  energy 
dissipation  in  a  structure. 

Four  main  methods  for  identifying  dominant  power  sources  and  power 
flow  paths  are  addressed  in  the  literature:  *  ;perimental  methods,  statistical 
energy  analysis,  the  finite  element  method,  an,d  a  finite  element  analogy  using 
heat  conduction  equations.  Experimental  solutions  are  the  most  common. 
The  authors  of  some  of  these  papers3-7  use  multiple  transducers  and  digital 
signal  processing  techniques  to  solve  various  power  flow  problems.  A 
common  method  is  the  calculation  of  cross  spectral  densities,  in  which  two 
accelerometers  are  placed  a  known  distance  apart  on  a  structure,  and  response 
spectra  are  generated  for  the  two  measurement  locations.  The  spectra  are 
statistically  analyzed,  and  power  flows  are  computed  over  some  range  of 
frequencies.  This  approach  is  similar  to  the  two-microphone  technique  used 
by  acousticians  to  solve  noise  propagation  problems  in  fluid  media.  Once  an 
experimental  apparatus  is  set  up,  the  analyst  may  easily  vary  applied  loads  and 
loading  frequencies.  Unfortunately,  the  added  weights  and  inertias  of  the 
transducers  attached  to  an  experimental  structure  may  produce  accuracy 
problems. 

Statistical  energy  analysis  (SEA)  is  a  computational  method  used  to 
solve  energy  flow  problems  in  the  high  frequency  domain.  A  definitive 
reference  on  SEA  is  the  text  by  R.L.  Lyon.8  In  the  SEA  method,  large 
structures  are  split  into  smaller  subsystems;  a  modal  density  is  estimated  for 
each  subsystem  so  the  number  of  modes  in  a  given  frequency  band  can  be 
determined;  dissipation  loss  factors,  which  relate  energy  stored  to  power 
dissipated,  are  estimated  for  the  subsystems;  and  coupling  loss  factors,  which 
relate  differences  of  modal  energy  of  subsystems  to  power  flow,  arc  assigned 
to  the  junctions  of  the  subsystems.  The  energy  distribution,  power  flows,  and 
power  dissipations  are  then  computed. 

SEA  is  a  good  way  of  solving  for  the  average  response  of  structures  at 
high  frequencies;  however,  all  spatial  variations  of  the  power  flow  field  in  t he 
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substructures  remain  unknown.  A  more  discrete  method  must  be  used  to 
identify  specific  power  tlow  paths  through  a  structure.  Finite  element  analysis 
(FEA)  may  be  used  for  this  purpose,  but  is  cost  effective  only  for  low  to  mid¬ 
range  frequencies,  since  higher  mode  shapes  are  more  complex,  wavelengths 
are  shorter,  and  denser  finite  element  meshes  are  required  to  model  a  problem 
correctly.  Mickol  and  Bernhard11  succesfully  used  FEA  to  identify  power  flow 
paths  in  simple  beam  and  plate  structures  excited  at  low  frequencies. 

Recently,  some  scientists  have  proposed  a  new  method  to  solve  for 
power  flows  in  the  middle  frequency  range.  Cuschieri,10  Net'ske  and  Sung,11 
and  Wohlever  and  Bernhard1  have  been  studying  this  new  approach:  a  finite 
element  analogy  in  which  powers  are  substituted  into  a  heat  conduction 
equation  and  the  power  flow  field  is  found  directly. 

In  this  paper  FEA  is  used  to  solve  the  power  flow  problem.  For  lower 
modes  the  method  is  accurate,  models  are  simple  to  build  and  modify  using 
modern  modeling  software,  and  analysis  results  are  viewed  easily  using  post¬ 
processing  graphics  packages.  Since  we  are  not  aware  of  any  commercial 
software  which  contains  a  power  flow  capability,  the  formulations  and 
computer  methods  are  developed  here  for  NASTRAN.1-  a  widely  used 
program  for  structural  dynamics.  The  FEA  studies  presented  in  other 
literature  consider  only  the  contribution  of  flexural  wave  motion  to  power  flow. 
Other  motion  types  of  power  flow,  such  as  axial  and  torsional,  are  ignored.  In 
this  study,  all  types  of  power  flow  are  considered. 

First,  general  methods  and  formulations  for  power  flows,  mechanical 
intensities,  power  dissipations,  input  powers,  and  output  powers  are  developed 
for  global  models,  beam  elements,  plate  elements,  solid  elements,  and  scalar 
elements.  The  required  NASTRAN  solution,  the  algorithm  of  the  power  flow 
processor,  and  the  use  of  I-deas  Supertab  (for  post-processing  and  graphics) 
are  outlined.  Two  test  models  are  analyzed  to  verify  the  methods:  a  simple- 
truss  and  a  beam-stiffened  cantilever  plate.  Finally,  conclusions  about  the 
method  are  drawn  from  results  of  the  test  case  analyses,  and  some  thoughts 
about  future  directions  for  work  are  discussed. 

THE  FINITE  ELEMENT  SOLUTION 


General  Methods 

A  typical  power  flow  cycle  is  shown  in  Fig.  1.  The  figure  shows  an 
arbitrary  structure  mounted  to  a  connecting  structure  by  a  spring  and  damper 
coupling.  A  dynamic  load  is  applied,  and  energy  flows  into  the  structure  at  the 
load  point.  The  input  power  then  flows  through  the  structure  along  multiple- 
flow  paths  denoted  by  arrows  whose  lengths  represent  power  flow  magnitudes. 
As  the  energy  flows  toward  the  mounting,  it  is  dissipated  by  material  damping 
and  sound  radiation  into  a  surrounding  medium,  and  the  flow  arrows  shorten. 
The  flow  and  dissipation  processes  continue  until  the  remaining  energy  exits 
the  structure  through  the  mounting  and  flows  into  the  connecting  structure. 
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Though  only  one  power  entry  point  and  one  exit  point  are  shown  in  this 
drawing,  multiple  loads  and  mountings  may  exist.  A  classic  text  which 
describes  the  flow  of  structure-borne  sound  is  the  book  by  Cremer.  Heckl,  and 
Ungar.13 

The  structural  dynamics  problem  may  be  solved  using  NASTRAN. 
The  structure  may  be  modeled  with  various  element  types;  mountings  are 
modeled  with  scalar  spring,  damping,  and  mass  elements.  Constraints  and 
loads  are  directly  applied.  The  steady-state  response  for  the  model  is  solved 
for  a  given  excitation  frequency,  and  the  power  flow  variables  are  calculated. 

Power  Flow  and  Mechanical  Intensity 

Power  is  defined  as  the  time  averaged  product  of  a  force  with  the  in- 
phase  component  of  velocity  in  the  direction  of  the  force.  For  time-harmonic 
analysis,  where  complex  numbers  are  used,  this  calculation  may  be  visualized 
as  taking  the  dot  product  of  the  force  and  velocity  phasors. 

Multiplying  one  complex  number  by  the  in-phase  part  of  another 
complex  number  is  the  same  operation  as  multiplying  the  first  number  by  the 
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complex  conjugate  of  the  other  number  and  taking  the  real  part  of  the  result. 
Therefore  a  general  formula  for  power  flow  in  a  structure  is 

Power  =  Real  [F  V  ],  (1) 

where 

F  =  force  vector  and 

"v*  =  complex  conjugate  of  velocity  vector. 

The  real  part  of  power  flow  is  called  the  active  power. 

Mechanical  intensity  is  power  flow  per  unit  area,  or  stress  (force  per 
unit  area)  multiplied  by  the  complex  conjugate  of  velocity.  Mechanical 
intensity  is  similar  to  acoustic  intensity  in  a  fluid  medium,  which  is  the  product 
of  pressure  and  the  in-phase  component  of  particle  velocity. 

Damping  and  Power  Dissipation 

Power  may  be  dissipated  in  several  ways:  by  material  damping,  by 
mountings  and  other  passive  and  active  devices,  and  by  radiation  as  sound. 
This  section  discusses  power  dissipation  due  to  damping.  For  now  only 
material  damping  is  considered  in  the  dissipation  process.  The  effects  of 
sound  radiation  will  be  considered  in  the  future. 

Power  dissipation  is  calculated  differently  from  power  flow  and  power 
input.  Since  power  dissipation  is  the  rate  of  energy  dissipation,  the  energy 
level  of  a  given  element  is  calculated  and  multiplied  by  its  damping  coefficient. 
Multiplying  the  energy  dissipation  by  the  angular  frequency  of  excitation  gives 
the  power  dissipated  in  that  element. 

The  effects  of  the  material  damping  coefficient  are  significant.  As  the 
damping  coefficient  is  increased,  the  power  dissipated  increases.  If  the 
damping  coefficient  is  zero,  no  scalar  damping  elements  are  applied  to  the 
structure,  and  no  sound  radiation  is  considered,  power  dissipation  will  be  zero 
and  no  power  flow  will  exist.  This  is  because,  with  no  damping,  forces  and 
velocities  will  be  exactly  90  degrees  out  of  phase,  and  the  in-phase  part  of 
velocity  with  force  is  zero.  Although  this  situation  is  physically  unrealistic,  it 
may  occur  in  a  finite  element  analysis. 

To  solve  for  power  dissipation,  energy  dissipation  must  first  be 
calculated.  The  energy  level  in  an  element  is  the  sum  of  the  element’s  kinetic 
energy  and  potential  energy.  Since  this  is  a  steady-state  problem,  and  the 
energy  is  a  time-averaged  quantity,  it  may  be  calculated  as  twice  the  kinetic 
energy: 

E  =  mw*,  (2) 

where 

E  =  energy, 

m  =  element  mass,  and 
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v  =  velocity. 


Power  dissipation  is  then  calculated  as 

Pd.ss  =  2rtrjE,  (3) 

where 

f  =  rotational  frequency,  and 
t]  =  material  damping  coefficient. 

The  rj E  term  is  the  energy  dissipation,  and  multiplying  by  the  angular 
frequency  gives  the  energy  dissipation  per  unit  time,  or  power  dissipation.  The 
result  will  be  a  real  number,  since  the  energy  calculation  multiplies  velocity  by 
its  complex  conjugate. 

Eq.  3  is  based  on  two  assumptions:  rj  is  small  (<  0.4),  and  f  is  a 
resonance  frequency.  The  calculation  of  power  dissipation  is  more 
complicated  for  off-resonance  frequencies.  For  a  one  degree  of  freedom 
(DOF)  system  the  equation  for  power  dissipation  is  multiplied  by  Jn/w,  where 
is  the  natural  frequency  and  w  is  the  excitation  frequency.  For  multi-DOF 
systems,  such  as  finite  element  models,  total  power  dissipation  for  a  given 
frequency  must  be  calculated  as  the  sum  jf  the  contributions  from  each  of  the 
model’s  mode  shapes.  This  requires  a  modal,  rather  than  direct  frequency 
response  analysis.  Although  a  modal  solution  would  generate  the  data 
necessary  for  accurate  power  dissipations,  error  may  be  introduced  to  the 
power  flow  computations.  Since  the  forces  and  velocities  of  a  modal  analysis 
are  calculated  as  the  sums  of  the  contributions  of  the  mode  shapes,  they  may 
be  inaccurate  if  an  insufficient  number  of  modes  are  used  for  the  summation. 
At  this  point,  power  flows  are  considered  the  more  important  of  the 
calculations  and  a  direct  frequency  response  analysis  is  performed.  Therefore, 
only  the  power  dissipations  calculated  at  or  near  resonance  frequencies  are 
accurate  using  Eq.  3. 

Since  the  calculation  of  power  dissipation  involves  the  element  mass, 
the  calculation  is  mesh-dependent.  As  mesh  density  increases,  element  power 
dissipations  decrease.  For  example,  if  a  beam  element  is  subdivided  into  two 
beam  elements,  the  original  power  dissipation  is  split  between  the  two  new 
beams.  Power  dissipation  can  be  made  independent  of  the  mesh  by  dividing 
by  element  volumes,  but  at  this  point  the  actual  power  dissipations  are 
calculated  because  of  their  importance  in  checking  power  balances  (sec  the 
Power  Balance  section  below). 

Also,  power  dissipation  is  directly  related  to  the  mode  shapes  of  an 
analysis,  so  areas  of  large  displacements  and  velocities  will  be  large  energy 
sinks,  and  nodes  (points  of  near  zero  displacement)  will  dissipate  almost  no 
power. 
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Power  Input 

Power  inputs  are  calculated  by  multiplying  input  forces  by  the  complex 
conjugates  of  the  resulting  velocities  at  the  load  points.  Total  input  power  is 
calculated  as 

Pin  =Real[VFiV’],  (4) 

i-t 

where 


i  =  load  point,  and 
n  =  number  of  loads. 

This  is  a  global  calculation  which  is  independent  of  element  type. 

Presently,  only  force  inputs  ere  considered.  Other  load  types  applied  to 
a  finite  element  model  such  as  displacements,  velocities,  and  accelerations, 
may  input  power  to  a  .tructure.  The  calculation  of  input  powers  for  these  load 
types  will  be  derived  in  a  subsequent  paper. 

Power  Output 

Power  output  is  the  power  that  leaves  the  system  through  its  mountings 
and  enters  the  connecting  system(s).  The  external  system  is  modeled  using 
spring,  damper,  and  mass  elements.  These  scalar  elements  must  be  connected 
to  additional  grid  points,  which  are  grounded.  The  forces  of  constraint  are 
combined  with  the  velocities  of  the  grid  attached  to  the  scalar  element  to 
calculate  power  output.  The  power  output  is  calculated  as 

Pout  =  Real  (SFjVj-J,  (5) 

j-i 

where 

j  =  grounded  grid,  and 
n  =  number  of  grounded  grids. 

Power  Balance 

The  previously  defined  terms  (power  input,  power  dissipation,  and 
power  output)  are  all  used  to  verify  a  power  balance  for  a  given  problem.  The 
power  balance  equation  is 

#elem 

Pin  =  Yj  Pdiss  +  Pout.  (6) 

k-1 

This  same  equation  is  used  in  SEA  theory.  Since  Pin,  PlH„,  and  Pdiss  are  all 
calculated  independently,  if  the  power  balance  equation  holds,  then  the  power 
flow  solution  is  correct  (assuming  the  original  finite  element  solution  is 
accurate).  This  power  equilibrium  equation  is  therefore  an  important  check 
on  the  power  formulations  and  calculations. 
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Element  Formulations 


Beam  Elements 

Much  of  the  literature  in  the  field  of  power  flow  is  devoted  to  beams. 
Noiseux14  and  Verheij15  described  methods  of  measuring  the  flexural  power 
flow  in  beams.  Other  authors,  such  as  Li,16  Wohlever  and  Bernhard.1  and 
Nefske  and  Sung,11  have  developed  power  flow  capabilities  for  beam  elements. 

All  the  methods  developed,  however,  consider  only  flexural  power  flow. 
Though  flexure  is  arguably  the  dominant  response  in  a  beam,  axial  and 
torsional  response  can  also  be  important.  For  this  reason,  all  possible 
components  of  power  flow  will  be  considered  here. 

Power  flow  methods  for  NASTRAN  BAR  element  types  are  derived  in 
this  section.  Either  lumped  mass  or  coupled  mass  solutions  may  be  used. 
Unfortunately,  torsional  inertia  for  the  BAR  element  type  is  not  calculated  by 
NASTRAN.  Concentrated  mass  elements  with  beam  torsional  inertias 
entered  as  masses  must  be  added  to  the  model  at  the  appropriate  DOF  to 
solve  for  accurate  torsional  power  flows. 

Power  Flow  and  Mechanical  Intensity.  A  diagram  of  the  BAR  element 
and  its  force  output  conventions  is  shown  in  Fig.  2,  where  plane  1  is  vertical 
and  plane  2  is  horizontal.  Since  a  beam  is  a  one-dimensional  element,  energy 
flows  in  only  one  direction:  in  the  local  x  direction,  or  along  the  length  of  the 
beam.  The  total  power  flow  for  a  beam  element  is 

Px  =  Real  [  -  (Fxv;+V1v;+V2v*+T^-M2^+M1^)j,  (7) 

where 

Fx  =  axial  force, 

V i  =  shear  force  in  y  direction, 

V2  =  shear  force  in  z  direction, 

T  =  torsion  about  x, 

M2  =  bending  moment  about  y, 

=  bending  moment  about  z, 

vj  =  translational  velocities  in  direction  i,  and 
=  rotational  velocities  about  axis  i. 

The  negative  sign  on  the  result  comes  from  force  and  displacement  direction 
conventions  for  the  element.  Negative  signs  appear  in  the  formulations  for  the 
plate  and  solid  elements  for  the  same  reason.  The  negative  sign  on  the  M: 
term  reflects  the  NASTRAN  force  output  convention.  In  Fig.  2,  M2  is  shown 
as  positive  in  the  opposite  sense  to  9y.  Therefore,  M20*  is  opposite  in  sign  to 
the  other  power  flow  components. 

Velocities  are  calculated  by  NASTRAN  for  each  grid  point,  and  beam 
forces  are  calculated  on  aa  •  ement  level.  This  difference  creates  a  problem, 
because  some  way  of  sob  '<  ‘or  a  power  flow  on  an  element  level  is  required. 
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Fig.  2.  The  BAR  Element 


The  solution  is  to  solve  for  a  power  flow  at  each  grid  point  and  then  calculate 
the  average  quantity  for  an  element. 

Shear,  axial,  and  torsional  forces  are  constant  through  the  element  and 
are  the  same  for  each  grid  point.  Bending  moments  are  calculated  at  each  end 
of  the  beam  element.  Velocities  are  solved  for  at  a  global  level,  and  a 
coordinate  system  transformation  must  be  performed  to  find  the  velocities  at 
the  element  level.  Since  the  grid  coordinates  and  beam  orientation  vector  are 
given,  the  velocity  transformation  is  straightforward.  After  power  Hows  at 
each  element  end  are  calculated,  they  are  averaged  to  give  an  element  power 
flow.  Power  flows  are  then  transformed  back  to  NASTRAN’s  basic 
coordinate  system. 

Two  important  observations  may  be  made  about  power  How  in 
NASTRAN  beam  elements.  Since  power  How  is  one-dimensional  in  beams,  it 
is  independent  of  mesh  variations.  Increasing  mesh  density  or  varying  the 
mesh  pattern  will  not  affect  greatly  the  power  flow  results  (assuming  the  mesh 
is  dense  enough  to  model  accurately  the  wave  shapes  of  the  solution).  Also, 
since  power  flow  is  dependent  on  element  force  quantities  (axial  and  shear 


forces,  torsion)  that  are  discrete  values  at  each  element,  the  power  How  and 
mechanical  intensity  quantities  are  not  continuous  across  beam  element 
boundaries  (from  element  to  element). 

Power  Dissipation.  The  energy  of  a  beam  element  includes  both 
translational  and  rotatory  terms,  and  is  calculated  as 

E  =  m(vxv*+vyv'+vzv*)  +  IxxMx+Iyy^y-rIzz<Mz»  (8) 

where 

m  =  element  mass, 

Iyy>Izz  =  niass  moments  of  inertia  about  cross  section, 

=  polar  mass  moment  of  inertia  about  beam  axis, 
v;  =  local  translational  velocities  in  direction  i,  and 
9\  =  local  rotational  velocities  about  axis  i. 

The  element  energy  is  multiplied  by  2Trj£un/oj,  as  in  Eq.  3,  to  yield  power 
dissipation.  As  mentioned  earlier,  the  un/u>  term  complicates  the  calculation 
so  that  only  power  dissipations  at  or  near  resonances  will  be  accurate  at  this 
time. 


The  rotational  inertial  energies  are  generally  small.  However,  if  the 
beam  lengths  are  long  with  respect  to  the  cross  section,  the  mass  moments  of 
inertia  become  important.  In  the  case  of  torsion,  where  the  only  large 
displacement  is  rotation  about  the  beam’s  axis,  the  polar  mass  moment  of 
inertia  term  dominates  the  energy  calculation. 

Plate  Elements 

Many  Naval  structures  require  the  use  of  plate  elements  to  accurately 
model  them.  Since  the  beam  element  formulation  included  all  components  of 
power  flow,  power  flow  capabilities  for  quadrilateral  plate  elements,  which 
consider  both  flexural  and  membrane  effects,  are  developed. 

The  literature  concerning  plate  elements  is  growing,  and  publications  by 
Mickol  and  Bernhard,9  Williams  et  al,17  Koshiroi  and  Tateishi,18  Noiseux,14 
Fahy  and  Pierri,19  and  Cuschieri,20  investigate  mechanical  intensities  and 
power  flows  through  plate  structures.  As  with  the  literature  for  beams, 
however,  most  approaches  consider  only  flexural  effects. 

Power  Flow  and  Mechanical  Intensity.  A  diagram  of  a  QUAD  element 
and  its  force  and  stress  output  conventions  is  shown  in  Fig.  3.  The 
quadrilateral  element  is  two-dimensional,  and  power  may  flow  in  the  local  x 
and  y  directions.  The  power  flow  in  the  x  direction  is  calculated  as 

Px  =  Real  [  —  (Vxv*— Mx^y-f-Mxytf*+Fxv*-t-Fxyv*)J;  (9) 

the  power  flow  in  the  y  direction  is 

Py  =  Real  [  —  (Vyv*+My(ix—  Mxy^y+Fyvy+Fyxv’)], 
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Fig.  3.  The  QUAD  Element 


where 

transverse  shear  forces, 
bending  moments, 
twisting  moment  about  x  and  y, 
membrane  forces, 
membrane  shear, 

local  translational  velocities  in  direction  i,  and 
local  rotational  velocities  about  axis  i. 

The  negative  signs  on  the  Mx0y  and  Mxy0y  terms  conform  to  the  NASTRAN 
force  output  convention.  These  bending  moments  are  opposite  in  sense  to 
their  corresponding  rotational  velocities.  All  forces,  both  flexural  and 
membrane,  are  per  unit  width.  Therefore  the  units  for  power  Hows  in  QUAD 
elements  are  power  per  unit  length. 

As  was  the  case  for  beam  elements,  grid  velocities  here  must  be 
transformed  to  the  local  element  coordinate  systems  to  be  used  in  the  power 


V,,  Vy  = 

MX  My  = 
M*y  = 
Fx,Fy  = 
Fxy,  Fyx  — 
V;  = 


flow  calculations.  After  the  calculations,  the  power  flow  vectors  are 
transformed  back  to  NASTRAN’s  basic  coordinate  system. 

Since  the  quantities  used  to  calculate  power  flows  for  plate  elements 
(forces  per  unit  width  and  velocities)  are  both  mesh  independent,  power  flow 
per  unit  length  is  also  mesh  independent.  To  solve  for  mechanical  intensities, 
power  flow  is  divided  by  plate  thickness,  giving  power  per  unit  area. 
Mechanical  intensities  for  QUADs  are  also  mesh  independent. 

Power  Dissipation.  The  energy  of  a  QUAD  element,  when  only  scalar 
mass  terms  are  considered,  is 

E  =  m(vxv’+vyv’+vzvz),  (10) 

where 

m  =  element  mass,  and 

v;  =  local  translational  velocities  in  direction  i. 

Element  energies  are  multiplied  by  2TrjIun/u  to  calculate  power  dissipations. 
Power  dissipation  terms  are  calculated  at  each  grid  point  and  averaged  to  solve 
for  the  element  dissipation.  This  calculation  is  mesh-dependent,  since  power 
dissipation  is  directly  related  to  element  mass.  The  same  difficulties  that 
applied  to  BAR  elements  due  to  the  un/w  term  apply  for  QUAD  elements. 

Solid  Elements 

Since  literature  in  the  power  flow  field  comes  largely  from  the 
experimental  sector,  solid  elements  are  generally  not  considered.  An 
experimentalist  cannot  place  a  measuring  device  inside  the  material  of  a  solid 
structure.  The  paper  by  Pavic,7  however,  describes  a  method  for  measuring 
structural  surface  intensity.  His  method,  which  involves  placing  transducers 
on  various  surfaces  of  a  machinery  system  to  measure  two-dimensional 
mechanical  intensities,  may  be  extended  to  three  dimensions.  Since  a  finite 
element  code  has  no  restrictions  on  making  "measurements”  internal  to  a 
structure,  mechanical  intensities  may  be  calculated  throughout  a  solid  model. 

Power  Flow  and  Mechanical  Intensity.  For  the  BAR  and  QUAD 
elements,  force  output  is  provided  by  NASTRAN.  For  solid  elements 
(isoparametric  IHEXi  element  types),  stress  output  is  solved  for  at  grid  points 
and  at  the  element  centroids.  Pavic7  uses  stresses  and  velocities  in  his 
formulation  of  structural  surface  intensities.  His  formulas  for  mechanical 
intensities,  extended  to  three  dimensions,  are 
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lx  —  Roal  [  (o,sv!t+rxyVy+rsav2)], 

Iy  =  Real  [  -  (CTyv'+r^v'+r^v*)], 

Iz  =  Real  [  -  (crzv’+r^Vs+r^v*)], 

where 

Ix,Iy,Iz  =  global  mechanical  intensities, 

(Tx,<7y,<7z  =  normal  stresses, 
rxy>ryz>rxz  = shear  stresses,  and 

v;  =  global  translational  velocities  in  direction  i. 

Since  element  stresses  are  given  in  the  basic  coordinate  system,  the  velocities 
do  not  have  to  be  transformed  to  element  coordinate  systems  as  they  were  for 
the  BAR  and  QUAD  element  types.  Calculations  are  made  at  each  grid  point 
and  averaged  to  calculate  the  element  mechanical  intensity. 

Since  the  quantities  used  to  calculate  mechanical  intensities  for  solid 
elements  (stresses  and  velocities)  are  both  mesh  independent,  the  mechanical 
intensities  are  also  mesh  independent. 

Power  Dissipation.  The  power  dissipation  calculations  are  similar  to 
those  for  the  QUAD  elements  (Eq.  10).  Only  mass  and  translational  velocity 
terms  are  considered.  Power  dissipation  terms  are  calculated  at  each  grid 
point  and  averaged.  Again,  since  the  element  mass  is  directly  related  to  the 
element  energy,  power  dissipation  is  mesh-dependent.  Also,  the  term  in 
the  calculation  is  a  problem  for  solids  as  well  as  for  beams  and  plates. 

Scalar  Elements 

Scalar  elements  may  be  used  to  simulate  mountings  and  structures 
connected  to  the  finite  element  model.  Scalar  springs,  dampers,  and  masses 
may  be  used  to  model  stiffness,  damping,  and  mass  effects.  These  elements 
may  be  important  for  certain  analyses,  such  as  those  for  a  structure  not  rigidly 
mounted.  In  certain  cases,  power  may  flow  out  of  a  structure  into  an  isolator, 
which  will  absorb  much  of  the  energy,  or  into  a  surrounding  medium.  The 
accurate  modeling  of  boundary  conditions  must  therefore  include  scalar 
element  types. 

Power  flows  and  power  dissipations  are  not  yet  computed  in  scalar 
elements,  but  the  presence  of  external  stiffnesses,  dampers,  and  masses  may 
significantly  affect  the  results  in  the  structural  element  types. 

Computer  Methods 

A  flow  chart  of  the  solution  process  is  shown  in  Fig.  4.  NASTRAN's 
Rigid  Format  8  (Direct  Frequency  Response)  is  used  to  solve  problems  for  any 
combination  of  load  cases  and  excitation  frequencies.  The  model  information 
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Fig.  4.  Power  Flow  Solution  Process 


and  problem  solution  output  are  written  on  a  UT1  file,  which  is  used  as  input 
to  the  McPOW  (Mechanical  POWer)  program.  After  the  power  computations 
have  been  made,  power  flows,  mechanical  intensities,  and  power  dissipations 
are  written  on  two  output  files.  One  file  contains  a  tabular  listing  of  the  power 
flow  results;  the  other  file  is  formatted  as  input  to  the  I-DEAS  Supertab'1 
post-processor,  which  is  used  to  interpret  visually  the  results. 

NASTRAN  Solution 

Before  Rigid  Format  8  is  run,  an  eigenvalue  extraction  (Rigid  Format  3) 
can  be  performed  on  the  model  to  determine  the  resonant  frequencies  and 
their  corresponding  mode  shapes.  Power  flows  can  then  be  measured  at 
response  peaks,  and  the  dominant  type  of  power  flow,  such  as  flexural  or 
axial,  can  be  predicted  by  examining  the  mode  shapes. 

Power  Flow  Algorithm 

The  program  McPOW  is  composed  of  four  main  sections:  the  model 
information  section,  the  NASTRAN  output  section,  the  power  flow 
calculation  section,  and  the  output  section. 
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The  first  section  reads  the  model  information  data  blocks  from  the  UT1 
file.  The  NASTRAN  output  section  reads  the  analysis  information  data 
blocks  and  assigns  forces  and  stresses  to  element  variables,  velocities  to  grid 
points,  and  input  loads  to  grid  points. 

The  power  flow  calculation  section  first  calculates  input  powers  using 
the  input  loads  and  corresponding  grid  velocities.  Next,  grid  velocities  are 
assigned  to  elements.  Power  flows,  mechanical  intensities,  and  power 
dissipations  are  then  calculated  using  element  forces,  stresses,  and  the 
velocities  of  the  element  grids. 

The  output  section  writes  power  flow  information  on  two  files.  The  first 
file  contains  a  tabular  listing  of  the  solution  variables  for  each  subcase  and 
frequency;  the  second  file  is  a  data  file  in  I-DEAS  Universal  file  format. 

Post-Processing 

The  user  may  analyze  the  power  flow  output  in  two  ways:  by  inspecting 
the  listed  output,  or  using  I-DEAS  Supertab’s  post-processor  to  draw  contour 
plots  and  arrow  plots.  Analyzing  the  tabular  output  is  a  good  way  to  check 
power  balances.  Power  input  is  equal  to  power  dissipated  plus  power  output. 
However,  to  visualize  the  entire  power  flow  solution  in  any  reasonably  complex 
geometry,  a  good  graphics  post-processor  is  required. 

Color  contour  plots  can  be  used  to  display  power  flow  magnitudes  and 
power  dissipations.  Power  flow,  however,  is  a  vector,  and  arrow  plots  are 
needed  to  display  the  direction  of  the  flow.  Other  authors,  such  as  Heekl,:; 
and  Koshiroi  and  Tateishi,18  have  used  arrow  plots  to  show  power  flows  in 
plate  structures.  An  alternative  unavailable  in  I-DEAS  Supertab  is  a 
combination  of  a  contour  and  arrow  plot,  which  can  illustrate  magnitude  and 
direction. 


TEST  CASES 

The  test  problems  described  here  illustrate  the  use  of  beam,  plate,  and 
scalar  elements.  Solid  elements  have  not  been  tested  yet. 

Simple  Truss 

Problem  Statement 

A  diagram  of  a  simple  truss  is  shown  in  Fig.  5.  The  factor  q  is  the 
assumed  loss  factor  of  the  material.  The  truss  members  are  constructed  of 
three  different  types  of  cross  sections.  The  model  was  attached  to  ground  at 
its  top  and  bottom  by  springs  in  all  six  DOF.  The  scalar  elements  simulated 
the  effects  of  fasteners  and  the  surrounding  structure(s).  An  end  force 
excitation  was  applied  in  all  six  DOF  over  a  range  of  frequencies.  The 
properties  of  the  two  W  type  sections  are  given  in  civil  engineering  handbooks. 
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Fig.  5.  Simple  Truss  Problem 


The  finite  element  model  consisted  of  74  BAR  elements,  with  each 
beam  section  having  a  different  mesh  density.  Section  1  was  modeled  with  40 
elements  of  0.25  m  length,  Section  2  consisted  of  10  elements  of  0.5  m  length, 
and  Section  3  was  made  up  of  24  elements  of  about  0.325  m  length.  For  the 
scalar  elements,  spring  constants  were  set  at  about  100  to  1000  times  the 
stiffness  of  the  members  at  the  appropriate  DOF;  the  loss  factor  for  the 
springs  was  set  at  10  times  the  material  damping  constant  or  0.2.  This  model 
is  a  good  general  test  of  the  power  flow  methods  previously  outlined,  since  it 
has  a  varying  mesh  density,  and  the  three  sections  have  different  beam 
properties. 

Results 

The  first  analysis  performed  on  the  model  was  an  eigenvalue  extraction 
(Rigid  Format  3).  Although  there  is  damping  in  the  model  and  the  modes  are 
actually  complex,  real  modes  may  be  calculated  to  estimate  the  resonances. 
The  first  50  modes  ranged  in  frequency  from  1.87  Hz  to  174  IIz.  The  three 
members  experienced  different  types  of  motion  in  each  mode  (i.e..  axial, 
flexural,  or  torsional),  showing  the  need  for  the  calculation  of  all  types  of 
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Power  (W) 


power  flows  in  beams. 

The  end  load  was  then  applied  for  frequencies  ranging  from  1  to  100 
Hz,  with  a  resolution  of  1  Hz.  The  plot  in  Fig.  6  shows  the  response  of 
section  three  at  100  Hz,  or  the  29th  mode  of  the  truss.  The  right  end  of  the 
plot  is  the  loading  point,  and  the  left  end  is  the  junction  with  Section  1  and  the 
mounting  at  the  bottom.  Both  power  flow  and  power  dissipation  are  plotted. 
Power  flow  decreases  as  it  propagates  along  the  beam  due  to  power 
dissipation.  Power  dissipation  oscillates  between  low  and  high  points, 
approximating  the  mode  shape  of  the  beam.  When  dissipation  is  large,  power 
flow  slopes  downward;  when  dissipation  is  small,  power  flow  remains  level. 

Powers  for  Trass  Saation  3,  f=100  Hz 


Fig.  6.  Power  Flows  and  Dissipations  for  a  Single  Frequency 

Although  the  type  of  plot  shown  in  Fig.  6  is  an  effective  method  of 
displaying  the  power  flow  response  for  a  specific  case,  spectra  plots  are 
required  to  illustrate  responses  over  the  entire  frequency  range.  The  set  of 
four  plots  in  Figs.  7  and  8  shows  power  flow  at  different  locations  on  each  of 
the  truss  sections.  Section  one  is  split  into  two  graphs,  one  is  for  the  top  half 
of  the  beam,  and  the  other  is  for  the  bottom  half. 

Power  flows  are  plotted  for  three  locations  along  each  member:  at  the 
beginning,  middle,  and  end  (d/L  =  0.0,  0.5,  and  1.0  respectively).  For  the  top 
and  bottom  halves  of  Section  1,  the  beginning  of  the  section  is  at  the  joint  with 
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Section  2;  tor  Sections  2  and  3.  the  beginning  of  the  member  is  at  the  load 
point. 

The  expected  response  is  a  relatively  uniform  lowering  of  each  curve  as 
power  flow  progresses  from  beginning  to  end  along  each  beam.  This  is  indeed 
the  case  for  some  frequencies.  Ho%vever.  at  some  joints,  such  as  the  junction 
of  Sections  l  and  2  and  the  junction  between  Sections  1,  3.  and  ground,  power 
flows  in  ways  that  are  less  intuitive.  As  a  result,  some  of  the  plots  "cross  over" 
each  other  and  power  tlow  increases  from  beginning  to  end. 

For  different  frequency  bands,  where  different  mode  shapes  are 
dominant,  the  nature  of  the  power  tlow  through  the  model  will  vary.  The  three 
power  How  diagrams  of  Fig.  9  show  some  of  the  wavs  in  which  power  may  How 
through  the  truss  model  in  this  analysis. 


Fig.  9.  Power  Flow  Diagrams  for  Truss  Problem 


Diagram  1  shows  power  entering  Sections  2  and  3  at  the  load  points  and 
(lowing  out  toward  Section  1.  At  the  junction  of  Sections  1  and  2.  the  power 
tlows  from  Section  2  into  the  upper  half  of  Section  l.  At  the  junction  o| 
Sections  1  and  3,  power  (lows  from  Section  3  into  the  scalar  elements 
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connected  to  ground  and  into  the  bottom  half  of  Section  1.  Power  then  Hows 
from  the  bottom  of  Section  1  to  the  top  of  Section  1,  then  up  to  the  scalar 
elements  at  the  top  and  out  of  the  model. 

Diagram  2  shows  a  similar  case,  but  with  two  differences.  At  the 
junction  of  Sections  1  and  2,  power  from  Section  2  Hows  into  the  top  and 
bottom  halves  of  Section  1;  at  the  junction  of  Sections  1  and  3,  power  Hows 
only  into  ground. 

In  Diagram  3,  power  is  input  only  into  Section  2,  and  Hows  into  the  top 
and  bottom  halves  of  Section  1.  The  power  in  the  bottom  half  of  Section  1 
flows  down  to  the  junction  of  Section  1,  Section  3,  and  ground,  where  some 
power  flows  out  of  the  model  and  some  flows  up  into  Section  3. 

The  cases  shown  in  Diagrams  1  and  2  are  the  most  common  as 
indicated  by  examination  of  the  printed  output.  Other  possibilities  exist  but  do 
not  occur  often  for  the  range  of  frequencies  analyzed.  The  type  of  power  How 
diagram  which  occurs  for  a  given  frequency  may  be  determined  by  examining 
the  plots  in  Figs.  7  and  8.  When  power  flow  increases  in  going  from  d/L  =  0.0 
to  d/L  =  1.0,  then  power  has  entered  the  beam  at  d/L  =  1.0.  When  power  How 
decreases  in  going  from  d/L  =  0.0  to  d/L  =  1.0,  then  power  has  entered  the 
beam  at  d/L  =  0.0. 

Response  peaks  in  the  graphs  shown  in  Figs.  7  and  8  correspond  to 
different  types  of  motion  in  each  section.  Some  peaks  represent  flexural 
motion,  some  are  due  to  axial  response,  and  some  are  torsional.  A  power 
flow  algorithm  which  considers  only  flexural  response  would  give  incorrect 
answers  to  this  problem. 

Power  balances  (Pjn  =  Pout+  V)  Pj;5S)  were  accurate  at  the  resonance 
frequencies  where  power  dissipations  could  be  calculated  accurately.  No 
power  balances  were  performed  for  off-resonance  frequencies  since  power 
dissipations  could  not  be  calculated  accurately  due  to  the  term  in  the 

calculations.  Power  flows  at  the  truss  joints  balanced  as  well.  A  calculation 
based  on  Kirchhoff’s  law  for  currents  can  be  made,  with  power  Hows  in  the 
BARs  connected  to  the  junctions  analagous  to  currents. 

Beam-Stiffened  Cantilever  Plate 

The  analysis  of  ribbed  stuctures  combines  the  power  flow  methods  for 
beams  and  plates.  Nilsson23  used  SEA  methods  to  predict  the  transmission  of 
structure-borne  sound  through  ribbed  plate  models.  Iiere,  FEA  is  used  to 
calculate  the  low  frequency  response  of  a  beam-stiffened  cantilever  plate. 

Problem  Statement 

A  diagram  of  the  model  is  shown  in  Fig.  10.  Like  t lie  truss  model,  the 
cantilever  plate  model  was  attached  to  ground  at  its  end  by  springs  in  all  six 
DOF.  The  scalar  elements  simulated  the  effects  of  fasteners  and  the 
surrounding  structure(s).  A  uniform  end  load  was  applied  in  the  axial. 


Fig-  10.  Beam-Stiffened  Cantilever  Plate  Problem 


transverse  shear,  and  bending  directions.  A  12  x  30  mesh  of  QUAD  elements 
was  used  to  model  the  plate  and  two  sets  of  30  BAR  elements  modeled  the 
stiffeners.  The  BAR  elements  were  offset  relative  to  the  plates.  For  the 
scalar  elements,  spring  constants  were  set  at  about  100  to  1000  times  the 
stiffness  of  the  members  at  the  appropriate  DOF;  the  spring  loss  factor  was 
set  at  10  times  the  material  damping  constant  or  0.2. 

This  model  illustrates  the  power  flow  capability  for  plate  elements  and 
helps  to  test  further  the  beam  element  formulation.  Also,  the  power  balance 
equation  is  checked  for  multiple  element  types  in  a  model;  the  total  power 
dissipation  in  the  beams  plus  the  total  power  dissipation  in  the  plates  must 
match  the  difference  of  power  input  and  power  output. 

Results 

An  eigenvalue  extraction  of  the  model  showed  the  first  20  natural 
frequencies  ranging  from  7.5  to  271  Hz.  Loads  were  applied  to  the  model  for 
a  frequency  range  of  5  to  250  Hz  with  a  resolution  of  5  Hz. 


Fig.  11  shows  a  plot  of  input  powers,  summed  over  the  loaded  nodes; 
and  output  powers,  summed  over  the  scalar  elements  attached  to  the  wall. 
The  input  power  is  approximately  equal  to  the  output  power  for  low 
frequencies  (less  than  25  Hz).  At  higher  frequencies  power  dissipation  in  the 
elements  increases  and  the  power  output  curve  falls  below  the  power  input 
curve. 


0.0  50.0  100.0  150.0  200.0  250.0 

Frequency  (Hz) 

Fig.  11.  Input  and  Output  Powers  for  Cantilever  Plate 

A  plot  of  power  flows  in  one  of  the  beam  stiffeners  over  the  frequency 
range  is  shown  in  Fig.  12.  Since  the  model  and  the  loading  function  are 
symmetric  about  the  center  of  the  plate,  power  flows  through  both  beam 
stiffeners  are  the  same.  Curves  are  plotted  for  three  locations  along  the  beam, 
with  d/L  =»  0.0  at  the  load  point.  The  curves  resemble  those  in  Fig.  11,  but 
with  lower  power  levels.  Also,  the  power  curve  at  the  load  points  is  generally 
lower  than  the  power  curves  at  locations  further  along  the  beam. 

To  explain  the  behavior  in  the  beam  stiffeners,  graphical  (contour  and 
vector)  plots  of  the  entire  domain  are  needed  to  supplement  the  spectrum 
plots.  A  contour  plot  of  the  mechanical  intensity  magnitudes  of  the  pUue  and 
beam  elements  is  shown  in  Fig.  13,  and  a  vector  plot  of  mechanical  intensities 
is  in  Fig.  14.  Mechanical  intensities  are  plotted  rather  than  power  flows  since 
mechanical  intensities  are  mesh  independent  for  both  element  types.  The 
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Frequency  (Hz) 

Fig.  12.  Power  Flews  for  Three  Locations  Along  Beam  Stiffener 

beam  elements  are  illustrated  as  plates  in  the  diagram  so  their  results  may  be 
visualized.  Figs.  13  and  14  show  how  power  flows  through  the  model  at  7.5  Hz 
(the  first  resonance  frequency,  a  cantilever  mode).  Power  flows  into  the  model 
through  shell  elements  next  to  the  stiffeners.,  but  is  then  channeled  down  the 
beams,  so  the  power  levels  at  beam  locations  further  along  the  plate  are  higher 
than  those  at  the  load  point  for  this  frequency.  About  2%  of  the  input  power 
is  dissipated,  with  the  rest  flowing  out  of  the  model. 

In  this  case,  most  of  the  power  is  flexural.  At  other  frequencies, 
however,  other  types  of  power  flow  become  dominant,  such  as  power  flow  due 
to  twisting  moments  and  in  plane  membrane  forces.  Also,  the  interchange  of 
energy  between  the  beam  stiffeners  and  the  plate  will  depend  on  the  mode 
shapes  that  are  dominant  for  a  given  frequency  range.  This  shows  that  the 
nature  of  a  power  flow  field  is  frequency  dependent.  In  general,  no  single 
trend  can  be  found  to  explain  power  flow  spectra  over  all  frequencies. 

The  plots  in  Figs.  13  and  14  show  the  power  flow  field  at  7.5  Hz,  and 
explains  how  power  flows  and  interchanges  from  plate  to  beam  stiffener. 
However,  other  plots  are  needed  to  explain  phenomena  at  other  frequencies. 
With  specialized  graphics  software  being  developed  by  Richard  Van  Escltine  of 
Code  1843  at  DTRC,  these  plots  may  be  animated  over  frequency,  allowing 
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the  analyst  to  observe  the  power  flowing  through  a  model  in  real  time.  Using 
such  a  tool,  the  influence  of  certain  modes  on  the  power  flow  field  can  be 
observed  easily,  and  the  nature  of  power  flows  through  a  given  model  may  be 
better  understood. 


POWER  FLOW  MAGNITUDES,  FREQ 
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Fig.  13.  Mechanical  Intensity  Magnitudes,  f=*7.55  Hz. 

Power  balances  are  accurate  for  resonance  frequencies,  with  the  total 
power  dissipations  of  the  beam  and  plate  elements  matching  the  differences 
between  input  and  output  powers.  These  results  show  that  both  element  types 
may  be  used  accurately  in  a  single  model. 


SUMMARY  AND  FUTURE  WORK 

A  general  capability  for  the  calculation  of  power  flow  variables  (power- 
flow,  mechanical  intensity,  power  dissipation,  power  input  and  power  output) 
has  been  developed  for  use  with  the  finite  element  code  NASTRAN.  BAR, 
QUAD2,  QUAD4,  HEXA2,  IIIEXi,  MASSi,  and  ELASi  element  types  are 
currently  supported.  Unlike  most  of  the  studies  described  in  other  literature, 
this  study  considers  all  types  of  power  flows  -  flexural,  axial,  and  torsional  -  in 
the  element  formulations. 
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Fig.  14.  Mechanical  Intensity  Directions  f-7.55  Hz. 

The  results  of  the  test  problems  indicate  the  method  developed  here  is  a 
valid  way  of  predicting  the  power  flow  response  of  a  dynamically  excited 
system  at  relatively  low  frequencies.  The  results  also  show  that  the  power  flow 
paths  in  complex  models  are  frequency  dependent.  As  different  modes  and 
types  of  power  flow  become  dominant,  the  nature  of  power  flow  fields  changes. 
To  better  understand  power  flow  fields,  graphical  animations  of  power  flow 
contour  and  vector  plots  over  frequency  can  be  used. 

The  use  of  FEA  to  calculate  power  flows  is  accurate  and  economical  for 
the  lower  modes  of  a  mechanical  system.  However,  the  power  flow  results  can 
be  no  better  than  the  NASTRAN  results  on  which  they  are  based.  Good 
modeling  techniques  and  an  understanding  of  the  wavelength  sizes  of  a 
problem  are  required.  The  shorter  the  wavelengths,  the  denser  the  required 
mesh  will  be. 

Planned  future  work  is  extensive.  A  list  of  topics  follows. 

1.  The  computational  methods  presented  here  will  be  compared  to  the 
experimental  methods  developed  by  David  Taylor  of  Code  274  of  DTRC. 
Both  beam  and  plate  models  will  be  analyzed  and  compared. 
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2.  Energy  is  not  only  dissipated  by  material  damping,  but  is  also  lost  due  to 
radiation  into  fluid.  Ross24  discusses  the  effects  of  entrained  mass  of 
fluid  and  air  on  the  sound  radiation  of  power.  Some  way  of 
implementing  the  effects  of  entrained  mass  will  be  investigated. 

3.  Once  power  flow  paths  and  power  sinks  are  identified,  some  way  of 
reducing  Pout  should  be  found.  Redman-White  et  al.25  and  Cuschieri26 
have  both  discussed  active  force  control  as  a  means  of  reducing  power 
flow.  Additional  forces  are  applied  at  calculated  locations  which  alter 
the  response  of  the  mechanical  system  so  that  power  flow  is  reduced. 
Some  power  sources  actually  act  as  power  sinks  in  this  case,  absorbing 
energy  from  the  system.  This  method  is  currently  being  investigated 
using  experimental  methods  by  John  Clemens  of  Code  274  of  DTRC. 
Any  computational  analyses  will  be  compared  to  the  experimental  ones. 
Another  method  of  reducing  power  flow,  which  has  been  used  for  many 
years,  is  the  application  of  dampers.  A  few  well  placed  dampers  can 
absorb  much  of  the  energy  in  a  mechanical  system.  A  relatively  new 
technology  that  could  be  applied  to  the  problem  is  shape  optimization.  A 
given  structure  could  be  optimized  for  minimum  sound  radiation,  output 
power,  and/or  weight. 
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APPENDIX:  NASTRAN  INPUT  DATA 

Several  data  blocks  must  be  written  by  NASTRAN  on  the  UT1  output 
file  for  the  McPOW  program.  The  following  ALTER  statements  are  put  in 
the  Executive  Control  Deck: 


s 

$  THE  FOLLOWING  STATEMENTS  CORRESPOND  TO  THE  1987,88 
S  VERSIONS  OF  COSMIC  NASTRAN,  RF8 
$ 

ALTER  23$  AFTER  THE  TA1  MODULE 

OUTPUT2  CASECC,EST,MPT,EQEXIN$ 

ALTER  135$  AFTER  THE  SDR2  MODULE 

OUTPUT2  OPPCl,OESCl,OEFCl,OUPVCl,OQPCl$ 

ENDALTER$ 


In  this  alter,  the  CASECC  data  block  contains  case  control 
information,  the  EST  data  block  holds  element  information,  the  MPT  data 
block  contains  material  properties,  and  the  EQEXIN  data  block  holds  grid  and 
SIL  (Scalar  Index  List)  information.  The  OPPC1  data  block  contains  the 
applied  forces,  the  OESC1  data  block  lists  the  element  stresses,  the  OEFC1 
data  block  holds  element  forces,  the  OUPVC1  data  block  contains  grid  point 
velocities,  and  the  OGPC1  data  block  contains  constraint  forces. 

To  ensure  that  all  the  required  data  are  in  the  data  blocks,  the  following 
output  requests  must  be  made  in  the  case  control  deck: 

FORCE=ALL 

STRESS=-ALL 

VELOCITY=ALL 

OLOAD=ALL 


Either  PHASE  or  REAL  output  data  may  be  specified. 
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